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Baryons in the large N limit of the massive two-dimensional 
Nambu— Jona-Lasinio model 
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Baryons in the massive Nambu-Jona-Lasinio model in 1+1 dimensions (the massive chiral Gross- 
Neveu model) are studied in the limit of an infinite number of flavors. The baryon mass is evaluated 
for a wide range of bare fermion masses and filling fractions, combining analytical asymptotic ex- 
pansions with a full numerical Hartree-Fock calculation. 
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I. INTRODUCTION 



Since the seminal paper by Gross and Neveu in 1974 
field theoretic models of fermions in 1+1 dimensions 
with (broken or unbroken) chiral symmetry have turned 
out to be quite useful. In view of a long list of appli- 
cations in particle, condensed matter and mathemati- 
cal physics, it is surprising that after more than three 
decades the basic solution of these models is still incom- 
plete, even in the tractable limit of an infinite number 
of flavors N. A few years ago it has been realized that 
the alleged phase diagram of Gross-Neveu (GN) mod- 
els in the^temperature (T) versus chemical potential (/x) 
plane [2, 3] was inconsistent with the known baryon spec- 
trum. In the standard GN model with scalar-scalar in- 
teraction, this problem has been solved in the meantime 
with the construction of a soliton crystal phase (see [1] 
and references therein) built out of Dashen-Hasslacher- 
Neveu baryons [5|. However there are still large gaps 
in our understanding of the model with continuous chi- 
ral symmetry, i.e., the Nambu-Jona-Lasinio model [6( in 
two dimensions (NJL2). In the chiral limit, it is known 
that solitonic baryons and soliton crystals do exist here 
as well, although their properties are quite distinct from 
those of the standard GN model In particular, in- 
duced fermion number plays an important role in the 
NJL2 model [8|. By contrast, in the massive NJL2 model, 
very little is known yet about either baryons or the phase 
diagram. We are aware of only two works addressing the 
issue of baryons at finite bare fermion mass. Salcedo et 
al. Q use an elegant variational ansatz to reduce the 
problem to the sine-Gordon kink, both in the 't Hooft 
model [lOJ and the NJL2 model, and point out the close 
relationship to the Skyrme picture in higher dimensions 
In Ref. [HJ the sine-Gordon equation was subse- 
quently identified as the first term of a systematic chiral 
expansion and higher order corrections were determined 
with the help of derivative expansion techniques. In the 
present work, we try to further improve this situation by 



constructing baryons in the (large N) NJL2 model for 
arbitrary bare fermion mass, using a combination of ana- 
lytical and numerical methods. Our ultimate goal is the 
phase diagram of the model which will be discussed else- 
where. Like in the GN model, we expect the baryon mass 
to play the role of the critical chemical potential where a 
soliton crystal phase sets in at zero temperature. Aside 
from their importance for the phase diagram, baryons 
in the massive NJL2 model are also of interest in their 
own right. By varying the bare fermion mass, we have 
the unique possibility to interpolate between the Skyrme 
picture characteristic for the chiral limit and the non- 
relativistic valence picture expected to hold in the heavy 
fermion limit. 

Let us now pose the baryon problem at large N in a 
more precise fashion. The Lagrangian of the (massive) 
NJL 2 model reads 

C = $(& - m )iP + 9 - + feVO 2 ] (1) 

where flavor indices are suppressed as usual, i.e., t/np = 
12^=1 ^kiftk etc. In the large N limit the Hartree-Fock 
(HF) approximation 13|, or, equivalently, the station- 
ary phase approximation to the functional integral [f|, 
become exact. The Dirac-HF equation 

(-75^ + ~/°S{x) + rfP(x)) i> a = E a il> a (2) 

for the single particle orbital a has to be solved simulta- 
neously with the self-consistency conditions 



m 



P = -Ng 2 J2i> a ij 5 ip a . 

a 

The HF energy 

OCC . - 



(3) 



(4) 
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comprises the sum over single particle energies of all oc- 
cupied states (including the Dirac sea) and the standard 
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correction for double counting of the potential energy. 
The vacuum problem is a special case of Eqs. (J2j{4j) where 
S = m (the physical fermion mass), P = 0, and all nega- 
tive energy states are filled. Self-consistency then yields 
the gap equation (in units where m = 1, and using the 
ultraviolet (UV) cutoff A/2 \v^) 



Ng 2 



7 + In A. 



(5) 



Here we have introduced the "confinement parameter" 7 

MM, 



7 = toq In A, 



(6) 



where the name (borrowed from condensed matter 
physics) refers to confinement of kink and antikink, not 
of the elementary fermions. Trading the bare parame- 
ters (g 2 , too) for the physical parameters (to = 1, 7) with 
the help of Eqs. (|5l6p is all that is needed to renormalizc 
the model and eliminate divergences in the baryon prob- 
lem. Baryons in the HF approach are characterized by 
x-dependent potentials S(x),P(x) and one (partially or 
fully) occupied extra level as compared to the vacuum. 
The observables of most interest to us are the baryon 
mass, defined as the difference in HF energies between 
baryon state and vacuum, the self-consistent potentials 
S, P and the fermion density. 

We shall use both analytical and numerical methods, 
depending on the parameters. The paper is organized 
accordingly. In Sec. HH we remind the reader of the 
derivative expansion and collect the results relevant for 
the vicinity of the chiral limit. In Sec. IHII we apply a 
recently developed effective field theory for the valence 
level [16[ to get complementary analytical insight into 
the non-relativistic regime. The full numerical solution 
of the HF problem including the Dirac sea is presented 
in Sec. IIVI Sec. [V] contains a short summary and our 
conclusions. 



II. DERIVATIVE EXPANSION 

In general the HF problem as defined through 
Eqs. (12IM[) is rather involved. This raises the question 
whether one can bypass its full solution, at least in some 
regions of parameter space. Near the chiral limit, the 
potentials S, P become very smooth. It is then possible 
to "integrate out" the fermions approximately, using the 
derivative expansion technique, see e.g. Refs. [13, EH- 
This method presupposes full occupation of each level 
with A" fermions and is inapplicable for partially occupied 
levels. It results in a purely bosonic effective field theory 
for the complex scalar field $ = S — \P which, in the 
large N limit, can be treated classically. Hence, one gets 
direct access to the HF potential without need to solve 
the Dirac-HF equation self-consistently. The leading or- 
der reproduces the sine-Gordon approach of Ref. [9j . In 
Refs. |12l [l9| this program was carried through to rather 



high order in the derivative expansion, yielding a sys- 
tematic expansion in the ratio of pion mass to physical 
fermion mass. The field $ was computed in polar coor- 
dinates, 



$= (1 + A)e 2ix , 



(7) 



with A, x expressed in terms of the pion mass and the 
spatial variable £ = m,i. For the present purpose, we 
express everything in terms of the confinement parameter 
7, using the (approximate) relationship 
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To this order in 7, the baryon mass is given by 



M B = 
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whereas the result for the potentials translates into 
2 2 



(8) 
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S = 1- 
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(10) 



In the vicinity of the chiral limit, fermion number is in- 
duced by the topologically non-trivial field $ and thus 
only resides in the negative energy states. For the fermion 
number divided by N ("baryon number"), the derivative 
expansion yields 



B 



Ax 



1 



[x(°°) - x(-oo) 



fill 



relating baryon number to chiral U(l) winding number 
of the field $, much like in the Skyrme model. One can- 
not conclude from this result that the induced fermion 
density is given by x'/ 71 " though. As discussed in Ref. [|[, 
the fermion density can be determined from the equation 
for the partially conserved axial current (PCAC), 



= 2m if)i"f b il). 



Taking the expectation value of Eq. ([T 



(12) 

in the baryon 



state, using the relationship jg = e^ u j v (valid in two 
dimensions) and integrating over x, one finds the fermion 
density per flavor 



p=-— [ dx'P(x'). 

n J-00 



(13) 



An expansion in 7 then yields in the present case 
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FIG. 1: Thick solid line: induced fermion density (O, dashed sion coefficient T (rising curve) for baryon HF potential JTO]) 
line: leading order term x'/ 71 "! thin line: difference, according versus momentum (7 = 0.2). 



to derivative expansion (7 = 0.2). 



Judging from the convergence properties of the series in 
7, we expect expressions (|8l TT4"l) to be quantitatively reli- 
able up to 7 ~ 0.2. 

Let us use these results to exhibit some differences be- 
tween baryons in massive GN and NJL 2 models, respec- 
tively. First of all, in the GN model, there is no induced, 
but only valence fermion density. The negative energy 
states cancel exactly when taking the difference between 
the fermion density of the baryon and the vacuum [13]. 
Turning to the NJL2 model, the situation is different. In 
Fig. [T] we illustrate the induced fermion density for the 
case 7 = 0.2 and compare it with the naive expectation 
y 7 / 71 "- The difference integrates to 0, in agreement with 
Eq. (fTTj) . Another pertinent observation is the following. 
In the GN model, the scalar potential S is reflectionless 
for any static solution [j, [2(|. This has turned out to 
be instrumental for the exact, analytical solvability of 
the baryon problem. We can now easily check whether 
the same is true for the NJL2 model, at least in the re- 
gion of validity of the derivative expansion. By solving 
the Dirac scattering problem with potential (fTU)) , we find 
that this potential is not transparent, see Fig. [2] for the 
case 7 = 0.2. The transmission coefficient rises from 
at low energies to 1 at high energies. Similar results are 
obtained at other values of 7. This is at variance with 
claims in the literature (2jJ (see however Ref. [22|]) and 
already indicates that the baryon problem in the NJL2 
model is more challenging than in the GN model. 



III. NO-SEA EFFECTIVE THEORY 

Recently, an approximation has been devised which 
yields analytical insight into the regime where the baryon 
features one positive energy valence level, while still hav- 
ing smooth HF potentials. The idea is to integrate out 
the negative energy fermions (i.e., the Dirac sea) and de- 



rive an effective theory for the (positive energy) valence 
fermions only. The HF equation for the baryon then re- 
duces to a single non-linear (and in general non-local) 
Dirac equation. This is a tremendous simplification as 
compared to the full relativistic HF problem, although 
somewhat more involved than the derivative expansion 
of Sec. lU where the fermions have been integrated out 
altogether. For the derivation of the effective Lagrangian 
in the massive NJL2 model, we refer to Ref. [lg]. Here 
we concentrate on the solution of the non-linear Dirac 
equation and the determination of the parameter region 
where the proposed truncation is valid. 

We proceed as follows. The effective no-sea Lagrangian 
consists of a massive free Dirac part and interaction terms 
descending from the original scalar and pseudoscalar 
four- fermion interactions. As far as terms involving only 
scalar interactions are concerned, we keep all contribu- 
tions in Eq. (39) of Ref. [l6[ . The terms involving pseu- 
doscalar interactions require one additional approxima- 
tion. Starting point is the result for the massive NJL2 
model, Eq. (76) of [Hj], replacing the scalar effective cou- 
pling constant g^ s — ir/N by the correspon ding expres- 
sion in the massive model (cf. Eq. (72) of [f6|). If we 
insist on an analytical solution as we do in this section, 
we are only able to handle the case in which the non- 
locality due to one-pion exchange is of short range. We 
therefore assume that the mass term m\ = 47 dominates 
the inverse pion propagator and expand in the remaining 
two terms as follows, 



1 



"ir JV(l+ 7 ) 



^75 %j) 



(15) 
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N 



ml Nm%(l + 7) 



The result is added to the GN model effective La- 
grangian. At this stage it is difficult to predict the precise 
range of applicability of the approximation due to the un- 
avoidable issue of self-consistency. We will return to this 
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question once we have solved the non-linear Dirac-HF equation. The final effective Lagrangian reads (m = 1) 



v\y )vt 2Ar(i + 7 ) vv< ^ 247V (1 + 7 ) 2 v ^ A ^ 6JV 2 (1 + 7 ) 3 VVVy 

llTT (1+ 67/11) , rn2 T.v,T ( v 7T 2 (1 + T /2) 7T 3 (1+7/4) 



PWXWO - , u W)(#) + ^7rrrr(#) 



1 



1440iV (I+7) 3 v ^/v^ 12JV 2 (1+ 7 ) 4V q N 3 (1 + 7 ) 5 

7T — 7T — — 7T^ I — — 

^(V5i 75 ^) 2 - ^^(□Vd75V0(#75V>) + 2jV 2 7 2 (1 + 7) (fe^) 2 (#)- (16) 



By construction, it is to be used in the positive energy 
sector only, since the effects of the negative energy states 
are already encoded in the Lagrangian. For a single 
baryon, the ensuing Euler-Lagrange equation is the non- 
linear Dirac equation for the (normalized) valence spinor, 



(17) 

Assuming the valence level to be occupied with filling 
fraction v = n/N (where n is the valence fermion num- 
ber) and using r\ = nv for ease of notation (all 7r's dis- 
appear), the scalar and pseudoscalar potentials are ex- 
pressed self-consistently as (' = d x ) 
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VSo 



r}4 
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(I + 7) 12(l + 7) 2 2(1 + 7 ) 3 
1177(1 + 67/ll)4 v 77 2 (1 + 7/2) 
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2?7 3 (l + 7/4) S 3 

3(l + 7) 5 

VPo _ VPo_ 

7 47 2 7 2 (1 + 7) 



6(1 + j) 4 



2s'^ + (stf 



27 2 (l+7) : 
V 2 PoSo 



(18) 



through scalar (sq) and pseudoscalar (po) valence level 
condensates 



So = ^oV>o, Po = ^0175^0- 



(19) 



It seems hopeless to solve the complicated non-linear 
Dirac equation (|17H19p in closed analytical form. How- 
ever, since the Lagrangian p^|) is an approximate one, it 
is sufficient to solve the equation perturbatively in some 
formal expansion parameter. We choose to expand in v 
and postpone till later the discussion of the validity of the 
truncation scheme. In order to solve the Dirac equation, 
we then find that we first have to solve a simple non- 
linear differential equation, followed by a sequence of in- 
homogeneous, linear differential equations. The lengthy 
analytical computations can be done conveniently with 
computer algebra (we used MAPLE), therefore we skip 
the details and record only the final results. The baryon 
mass (computed from the classical field energy) has the 



following power series expansion in r\ = ttv, 



Mb 
N 



(7 2 - 157 - 8)t7 5 



' 24(l + 7) 2 19207(1 + 7 ) 5 
( 7 4 - 302 7 3 + 265 7 2 + 376 7 + 88)r) 7 
3225607 2 (1 + 7 ) 8 ' 



(20) 



The energy E$ of the valence level is closely related to 
Mb- According to standard HF theory, the single particle 
energy can be interpreted as removal energy of a fermion, 
or equivalcntly 



En 



8Mb 
dn 



(21) 



This is indeed what we find in the calculation, so that 
there is no need to spell out Eq. Next we turn to the 
HF potentials (TS")) derived from the self-consistent va- 
lence level spinor ip . For the scalar and pseudoscalar 
potentials, we obtain 



s 22 2 



s 44 



cosh £ 

S62 



V cosh 2 £ cosh £ 



S64 



S66 



cosh 2 £ cosh 4 £ cosh 6 £ 



P = 



P33 3 



cosh £ 



P53 



+ 



P55 



cosh 3 £ cosh 5 £ 



V 



(22) 
sinh £ 



with 7 -dependent coefficients s nm ,p nm relegated to the 
appendix in order to keep the paper readable (the sub- 
scripts n, m denote the powers of 77 and 1/ cosh£, respec- 
tively). The spatial variable £ = yx in Eq. (|22|) involves 
a scale factor 



V = 



V 



( 7 2 - 37 - 2)if 



2(1 + 7) 48 7 (l + 7 ) 4 
(3 7 4 - 226 7 3 - 65 7 2 + 687 + 24)r? 5 
11520 7 2 (1 + 7 ) 7 ■ 



(23) 



One verifies that it is related to the single particle energy 
Eq via 



£ = Vi-y 2 - 



(24) 
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This last relation can be understood in physics terms as 
follows. Asymptotically, the valence wave function falls 
off exponentially with the k value (imaginary momen- 
tum) of the bound state. Due to self-consistency, the 
same parameter will govern the asymptotic exponential 
decay of the potentials in a no-sea HF calculation. This 
is exactly what Eq. (|24|) . together with the shape of the 
potentials S and P, guarantees. 

We note in passing that the asymptotic expansions for 
the scalar potential S at small 7, Eq. (jTUJ) , and at large 
7 or small u, Eq. (|22[) . have the same functional form. 
This is not true for the pseudoscalar potential P where 
Eq. (fT0|) involves even, Eq. (|22p odd powers of l/cosh£. 

We have tested the above results in two different ways. 
If we switch off the pseudoscalar coupling, we can carry 
out the same calculation for the massive GN model and 
reproduce the well-known exact results to the expected 
accuracy. A more specific test of the NJL2 model calcu- 
lation is provided by the divergence of the axial current 
implying Q 

d x p=-^P. (25) 

Here p denotes the fermion density per flavor consisting of 
the valence density and the induced part from the Dirac 
sea (n — uN), 

p = vipliPo - ^-d x P. (26) 

Eqs. (|25l26p yield the following non-trivial identity relat- 
ing the valence fermion density to the pseudoscalar HF 
potential, 

vdM^) = - 2 - 1 P + ^-dlP. (27) 

In our calculation it is violated at order v 7 . For the sake 
of completeness, let us quote the separate results for the 
valence fermion density p VB \ — vipoipo, 

^ _ V22 2 , ( ^42 «44 \ 4 

cosh 2 £ V cosh 2 £ cosh 4 £ / 

+ ( Ve2 _|_ 1,64 + V ™ \ ^ 6 / 2g % 

V cosh 2 £ cosh 4 £ cosh 6 £ / ' 
and for the induced fermion density pi n d = — j-c^P, 

/ U2 «44 \ 4 
TTpind = -J- + "TT V 

\cosh £ cosh £/ 

C J«_ + + _ie^_V e (29) 

Vcosh 2 £ cosh 4 £ cosh 6 £/ 

where the coefficients can again be found in the appendix. 

Let us now discuss the range of validity of our trunca- 
tion. The most obvious candidate is the regime 

7 > 1 (regime I) (30) 



and arbitrary v. In this case we can read off from our 
results that 'tjnp ~ I/7, ?/>i75'0 ~ V7 2 an d D ~ 1/7 2 - 
Inspection of the coefficients in C c s then shows that we 
have kept all terms up to order I/7 8 in an asymptotic 
expansion for large 7. Since the identity (|2T|) is violated 
at order I/7 7 we can trust the results (expanded in 1/7) 
up to order I/7 6 . By way of example, we note for future 
reference the asymptotic behavior of Mb for 7 — > 00 and 
full occupation of the valence level (y = 1) complemen- 
tary to Eq. ©, 

M B tt 2 ^ _ 7r 2 (240 - 7r 2 ) 

N 24 7 2 I2 7 3 19207 4 1 ' 

7T 2 (16-7T 2 ) 7T 2 (67200- 137767T 2 + TT 4 ) 

+ 967 s 3225607 6 ' 

Another case of interest is weak occupation of the va- 
lence level, i.e., taking the formal expansion in v liter- 
ally. Here, iptp ~ v 2 , ipvysip ~ ^ 3 ,D ~ v 2 and all terms 
up to order v s are kept in £ c ff. A necessary condition 
for this to be valid is v <C 1. Due to the expansion of 
the propagator in Eq. (| 15[) . inverse powers of 7 appear 
in the coefficients, and we have to require in addition 
v 2 <C 7. In this regime, all terms are then consistently 
kept and, according to the test based on Eq. |27|) . the 
results are trustworthy up to order i/ 6 . This 2nd regime 
of applicability may be summarized as 

v -C min(l, ^7) (regime II). (32) 

If the condition v <C y/j does not hold, the non-locality 
of the effective action becomes long range and our ap- 
proximation breaks down. 

We wind up this section with a summary of the analyt- 
ical results obtained in Sees. ITT1 and IIIII From a practical 
point of view, it is fairly easy to judge the quality of 
an effective theory by comparing successive orders in the 
expansion. This is illustrated in Fig. [3] where we plot 
the baryon mass (per flavor) Mb/N as computed from 
Eqs. ([5]) and (|2H|) versus v and 7. We only display the 
result if the highest order correction kept contributes less 
than 0.002 of the total result (other numerical values for 
the tolerance would give qualitatively similar plots). This 
shows clearly that the derivative expansion (curve start- 
ing at v = 1,7 = 0) and the no-sea effective theory (tilted 
surface) are valid in disjoint regions of the diagram, sepa- 
rated by a gap. Needless to say, the no-sea calculation is 
also valid at higher values of 7 not shown here. The struc- 
ture of the baryons is very different, depending on which 
asymptotic expansion we look at. Baryons described by 
the derivative expansion have no filled positive energy 
level, so that baryon number is entirely due to fermion 
number induced in the Dirac sea. If on the other hand 
the no-sea effective theory applies, there has to be a posi- 
tive energy valence level, and the valence fermion density 
dominates over the induced one, cf. the powers of 77 in 
Eqs. (|28l29p . Our formulae also show that in the limit 
of large 7 for fixed ^, or in the limit of small v for fixed 
7, the pseudoscalar terms are suppressed and the NJL2 
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FIG. 3: Analytical results for the baryon mass as a func- 
tion of filling fraction v and confinement parameter 7. Shown 
are only results which are stable at the 0.2% level if one re- 
moves the highest order correction term. Curve starting at 
7 = 0, v = 1: derivative expansion. Tilted surface: No-sea 
effective theory. 



results converge to those of the GN model. In fact, a non- 
relativistic approximation to the Dirac equation would be 
perfectly adequate in this regime. The region where the 
uppermost filled level reaches the middle of the gap and 
crosses zero energy lies in the white part of Fig. [3] and has 
to await a numerical treatment, to which we now turn. 

IV. NUMERICAL DIRAC-HARTREE-FOCK 
CALCULATION 

A. Hamiltonian in discretized momentum space 

We first set up the Dirac Hamiltonian in the form of 
a matrix which can readily be diagonalized numerically. 
Consider the Hamiltonian 

H = To -^-+ 1 °S(x)+i 1 1 P(x) (33) 
1 ox 

with given scalar and pseudoscalar potentials. We choose 
the representation 

7° = err, 7 1 = -icr 2 , 75 = I ! 1 = 0-3 (34) 

of the 7 matrices leading to a real symmetric Hamiltonian 
matrix below. In order to discretize the spectrum of H, 
we work in a finite interval of length L and impose anti- 
periodic boundary conditions for fermions, resulting in 
the discrete momenta 

We decompose the potentials into Fourier series, 

S(x) = Sie i2 * ex / L , P(x) = iJ2Pic i27:ex/L - (36) 

1 i 



Guided by the analytical results at small and large 7 
of the preceding sections, we assume that parity is not 
spontaneously broken in the baryon. A proper choice of 
the origin in coordinate space then leads to the follow- 
ing symmetry relations for the scalar and pseudoscalar 
potentials, 

S(x) = S(-x), P(x) = -P(-x). (37) 

Together with the reality conditions for S and P, we then 
find that the coefficients Si, Pi are real and satisfy 

S- t = Si, P- e = -P e , P o = 0. (38) 

The zero mode So of S(x) acts like a mass term and will 
be denoted by to (not to be confused with the physical 
fermion mass in the vacuum set equal to 1 throughout 
this paper). We decompose the Hamiltonian as follows, 

1 d 

H = H + V, H = 7 5 -_+7 m , 

1 ox 

V = j°(S(x) -to) +ij 1 P(x), (39) 

and diagonalize the free, massive Hamiltonian Hq first, 

H \q,n) = r?\An + m 2 \r],n}. (40) 

Here, r\ = ±1 is the sign of the energy and k n the discrete 
momentum. The free eigenspinors are given by 

with the shorthand notation k = k n ,E = \/k 2 + to 2 . 
Evaluating matrix elements of V in this basis, we find 

tf,ri\V\T,,n) = — L= (A^St+A^Pt) 

(42) 

(k' = k n , , E' = ^{k'Y + m 2 ) with 

4 ±) = r)'y/E' + r]'k'y/E ~r 1 k± ij^E + rjk^E' - rj'k'. 

(43) 

As Hq is diagonal, we can now easily construct the full 
matrix (77', n'\H |to n). We truncate the basis at some 
momentum index N, keep positive and negative energy 
states with labels —N — 1, ...,N and diagonalize the re- 
sulting (47V + 4) x (4/V + 4)-dimensional, real symmetric 
matrix H numerically. 

B. Perturbation theory deep down in the Dirac sea 

It is possible and in fact necessary to treat the negative 
energy states deep down in the Dirac sea perturbatively. 
First, this restricts the Hamiltonian matrix to be diago- 
nalized numerically to manageable size. Secondly, there 
is a logarithmic UV divergence in the sum over negative 
energy states which gets cancelled by a similar divergence 
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in the double counting correction. This cancellation has 
to be accomplished analytically before we can hope to 
reliably extract the finite part. We therefore compute 
the eigenvalues of H in 2nd order perturbation theory in 
V, in the large volume limit L — ► oo. For the negative 
energy continuum states, we find 



+ k 2 P 2 - 2k f E{k)SfP e 
(k 2 - k 2 )E(k) 

. (44) 

with E(k) = \Jk 2 + to 2 and kg = tt£/L. This result is 
valid for k 3> 7r£ max /L, where the inverse denominators 
don't blow up. By comparing the perturbative eigen- 
values with the ones from diagonalization, we can check 
whether we have taken into account a sufficient number 
of basis states and identify any eigenvalues affected by 
the truncation. The result is an optimal energy E m in, 
the energy of the lowest single particle level in the Dirac 
sea computed by matrix diagonalization. 



contribution AE2 to the baryon mass, 
AE> = _. )Z ; 'dfc'g E(M) 2 Sf + kjP 2 - 2k i E(k)S i P t 



femin 



(k 2 - k 2 )E{k) 



(48) 

and can be computed analytically (fc min is determined by 
Em_i-a introduced above). The first two terms in Eq. |[4"T|). 
after vacuum subtraction, give rise to an elementary, log- 
arithmically divergent integral which has to be regular- 
ized with the same cutoff as the gap equation and will be 
denoted by AE 3 , 



AE, = -21 / ^( E (k) + J2 



Sj + Pf 

? 2.^ + ^^ 



Vfc 2 + i) 



(49) 



— (m 2 -l + 2£(^ + P|))lnA + A J B 3 



C. Computation of the baryon mass 

The baryon mass has to be evaluated relative to the 
vacuum, i.e., as a difference of the HF energies ^ of the 
baryon and the vacuum. In the NJL2 model (like in the 
standard GN model, for that matter) the self-consistent 
potentials in the vacuum are S(x) — \,P(x) = in ap- 
propriate units. The gap equation ([5]) serves to eliminate 
the bare coupling constant from the double counting cor- 
rection in Eq. For the present discussion it is helpful 
to split the calculation of the baryon mass into four dis- 
tinct pieces 



M B /N = AE 1 + AE 2 + AE3 + AE 4 



(45) 



defined as follows. We first evaluate the sum over single 
particle energies for occupied states resulting from matrix 
diagonalization (subtracting the vacuum), 



AE 1 =J2(E a - ED . 



(46) 



Eroin has been introduced before, £max denotes the en- 
ergy of the highest occupied level in the baryon. The 
vacuum single particle energies E™ c are computed in a 
finite interval (length L) from the free massive theory 
with m = 1. Next we decompose the perturbative single 
particle energies (|4"4")l according to their UV behavior, 



E pcrt = -E(k) 



E 

1=1 



pi 



E 

t=i 



E{k) 
E(k e ) 2 S 2 + k\Pj 



(47) 



2k e E(k)S e P i 



{k 2 - k 2 )E{k) 



The last term on the right-hand side behaves asymptoti- 
cally like 1/fc 2 and can be integrated over momenta with- 
out cutoff. The resulting convergent integral defines the 



AE3 is the finite part which can again be computed ana- 
lytically. Finally, AE4 is the (vacuum subtracted) double 
counting correction, 



AE 4 = — m 2 
2vr I 



1 + 2J2(S 2 + P?) (7 + In A) 



(to — lj. 

7T 



(50) 



In the sum AE3 + AE4 the In A terms are cancelled ex- 
actly and a finite result for M~b/N, Eq. ([45|) . is obtained. 



D. Finding the self-consistent potential 

In non-relativistic HF calculations, one usually deter- 
mines the self-consistent potential iteratively. Starting 
from some guess, one solves the HF equation, computes 
the new HF potential and repeats this procedure until it 
has converged. This method was also employed in Ref. Q 
for the 't Hooft model, using a lattice discretization in co- 
ordinate space. We do not have this option here. The 
reason is the fact that we work in the continuum and have 
already sent the UV cutoff to 00 and the bare coupling 
constant to 0. As is clear from the self-consistency rela- 
tions ([3]), in this case the HF potential is determined by 
the far UV region where lowest order perturbation theory 
holds. Therefore we always get back the potential which 
we put in, irrespective of what we choose. If we would 
keep the cutoff finite like on the lattice, we could in prin- 
ciple use an iteration, but the convergence would slow 
down with increasing cutoff. A way out of this problem 
was found in Ref. [23| for the GN model where the HF 
energy of the baryon was simply minimized with respect 
to the potential S(x). Technically, the minimization was 
done with respect to the Fourier coefficients St of S{x). 
In the NJL2 case there are twice as many parameters due 



to the pseudo-scalar potential, so that this method would 
be rather cumbersome. Moreover, owing to the analyt- 
ical calculations in Sees, [ill and we already have a 
fairly good qualitative idea of how the potentials should 
look like. Hence the following strategy has turned out 
to be more economic. We write down a parametriza- 
tion of S and P which is both flexible and well suited to 
the known limiting cases. We then vary with respect to 
the parameters. Since we are only interested in a lim- 
ited precision, a relatively small number of parameters 
is sufficient, provided the parametrization is judiciously 
chosen. The convergence can be checked by comparing 
runs with different numbers of parameters. A natural 
ansatz for trial functions in the present problem is given 
by 

(£ = yx). This is clearly compatible with the no-sea 
effective theory, see Eq. (f22j) . The derivative expansion 
would suggest even powers of l/cosh£ in P(x) at low 
7 rather than odd powers, see Eq. (fTU|) , but this makes 
almost no difference in practice, once we keep a sufficient 
number of terms. The coefficients Ck,dk and the scale 
parameter y serve as variational parameters. The baryon 
mass computed as explained above is a function of these 
parameters and is minimized in the (2K + l)-dimensional 
parameter space via some standard algorithm. 

E. Numerical results 

All the results shown below were obtained by keeping 
the first three terms in each sum, Eq. (|5ip . and varying 
with respect to 6 coefficients (cfc, dk) and the scale param- 
eter y. We found that a matrix dimension 804 x 804 (cor- 
responding to N = 200 in Sec. II V A]) and box size L = 20 
(in units where m = 1) were adequate. The minimum 
in the 7-dimensional space was found with a standard 
conjugate gradient method involving 70 iteration steps. 
We first tested the whole procedure for several values of 
v, 7 where the analytical approaches are expected to be 
reliable. For example, at 7 = 0.2, v — 1 the derivative 
expansion gives Mb = 0.52227 and the numerical HF cal- 
culation Mb = 0.52248. This is satisfactory, especially 
since the numerical value is the result of subtracting large 
numbers (of the order of 10 000) due to the Dirac sea, cf. 
Eq. (|4"rJ)l . The scalar and pseudoscalar HF potentials are 
indistinguishable on a plot if one compares the deriva- 
tive expansion with the full HF calculation, therefore we 
don't show them here. At 7 = 5.0, v = 1 the no-sea ef- 
fective theory yields Mb = 0.98847 as compared to the 
numerical HF calculation Mb = 0.98850, again with ex- 
cellent agreement of the self-consistent potentials. We 
also tested the HF calculation against the no-sea effec- 
tive theory at small filling. For 7 = 0.2, v = 0.2 the 
analytical value M-q = 0.197314 compares well with the 
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12 3 4 

7 

FIG. 4: Circles: numerical results for baryon mass at v — 1 
versus 7. Curves: analytical asymptotic predictions (deriva- 
tive expansion to the left, no-sea effective theory to the right). 



In 




(15 ' i ' 1.5 2 2^5 

7 



FIG. 5: Comparison of baryon mass for large -/V GN (up- 
per curve) and NJL2 (lower curve) models, at v = 1. These 
curves can be regarded equally well as zero temperature phase 
boundaries in the (fx, 7)-plane. 



numerical result Mb = 0.197337. These examples give 
strong support to the numerical method as well as to the 
(independent) analytical calculations. 

Regarding the phase diagram, the question of most 
immediate interest is the dependence of Mb on 7 for full 
occupation v = 1. This is expected to yield a critical (2nd 
order) curve in the phase diagram at T = 0. In Fig. [4] we 
show the numerical points together with the asymptotic 
expansions, using both derivative expansion and no-sea 
effective theory. In Fig. [5] we compare this newly calcu- 
lated baryon mass in the massive NJL2 model with the 
known one from the massive GN model, choosing a some- 
what smaller range of 7 as compared to Fig.[5]to highlight 
the differences. The two curves shown reflect directly the 
critical lines in the (fi, 7)-plane at T = of the respec- 
tive phase diagrams. As a by-product of the calculation 
of baryon masses we also get the single particle energy 
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FIG. 6: Energies of valence levels for GN model (dashed 
curves) and NJL2 model (solid curves) versus 7, illustrat- 
ing how the NJL2 baryon interpolates between a Skyrme-type 
baryon with induced fermion number and a conventional, non- 
relativistic, valence-type baryon. 



of the valence levels, i.e., discrete levels inside the mass 
gap. It is instructive to compare these energies between 
the GN and NJL2 models as well. Fig.|5]shows the pair of 
(charge conjugation) symmetric states in the GN model 
which come together at zero energy at 7 = (this point 
corresponds to a kink and antikink at infinite separa- 
tion) . The baryons in the N JL2 model also feature a pair 
of discrete states which approach those of the GN model 
at large 7. At small 7, they converge to the lower edge of 
the mass gap, the upper level crossing zero near 70 = 0.3. 
Both of these valence levels are fully occupied. Below 
70 one would talk about induced fermion number (the 
number of negative energy levels changes as compared to 
the vacuum), whereas above 70 one has ordinary valence 
fermions like in the GN model. Note however that noth- 
ing discontinuous happens at the "spectral flow" point 
7o- This picture nicely illustrates the transition from the 
massless baryon at 7 = (with baryon number given by 
winding number of the pion field, like in the Skyrmion 
case) to the weakly bound non-relativistic baryon emerg- 
ing at large 7 in either GN or NJL2 models. The latter 
one is apparently insensitive to the original type of chiral 
symmetry, be it U(l) or Z2. 

The HF potentials belonging to these calculations are 
shown in Figs. ITIfTOl Fig. [7] is a 3d-plot containing the 
full information how S and P depend on x for 20 differ- 
ent values of 7. The 3 different projections of this plot 
onto the coordinate planes show the x-dependence of S 
(Fig. [8]) and P (Fig. [9j) in more detail, as well as paramet- 
ric contour plots in the (5, P)-plane (Fig. [TO]). We can 
infer from this last figure that the chiral winding num- 
ber jumps from 1 to slightly below 7 = 0.7, where the 
contour plot hits the origin (close to the 5th curve from 
the outside). Note that this does not coincide with the 
value 70 = 0.3 where the upper valence level crosses zero, 
so that induced fermion number is not directly related to 
winding number of $ = S — iP for larger bare fermion 
masses. This is in agreement with the general observa- 




FIG. 7: HF potentials S, P versus x in a 3d plot. The curves 
correspond to 7 = 0.3, 1.0 in steps of 0.1 and 7 = 1.2, 3.4 
in steps of 0.2. The outermost curve has the smallest value 
of 7. The curves at large 7 accumulate around S = 1, P = 
and are not well resolved, see also Figs. [8 U10I 
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FIG. 8: Projection of Fig. [7] onto the (x, S)-plane, showing 
how the shape of the scalar potential evolves with 7. 



tion that the integer part of induced fermion number is 
not topological, being sensitive to spectral flow [24I ]. 

Although it is presumably not relevant for the phase 
diagram, we now turn to baryons with partially occupied 
valence level. The HF calculation can be done exactly 
as before, except that the single particle energy of the 
upper valence level is weighted by the occupation fraction 
v when calculating the baryon mass. We have carried out 
detailed studies of the ^-dependence for the three values 
7 = 0.2,0.6,1.0 and compared the HF results with the 
no-sea effective theory expected to describe correctly the 
v — ► behavior. The results for the baryon mass are 
shown in Figs. 1X111X31 and again underline the value of the 
much simpler analytical results. 



10 





FIG. 9: Analogous to Fig. [8] but for the (x, P)-plane and 
pseudoscalar potential. 



FIG. 12: Same as Fig. [TT] but for 7 = 0.6. 
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FIG. 10: Projection of Fig. E]onto the (S, P)-plane. 
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FIG. 11: Baryon mass in the NJL2 model versus filling frac- 
tion v of upper valence level, for 7 = 0.2. Circles: numerical 
HF calculation, curve: analytical prediction based on the no- 
sea effective theory. 




FIG. 13: Same as Fig. EH but for 7 = 1.0. 




FIG. 14: HF potentials S, P versus x in a 3d plot. The 17 
curves correspond to v = 0.2, 1.0 in steps of 0.05. The 
outermost curve has the largest value of v. 
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FIG. 15: Baryon mass as a function of filling fraction v and 
confinement parameter 7. The numerical HF results (circles) 
have been combined with the analytical results of Fig. [3] com- 
pleting the picture of what is known by now. 



The self-consistent potentials as a function of v are 
illustrated for the case 7 = 0.2 in Fig. [14] (the other cases 
are qualitatively similar). This 3d-plot corresponds to 
Fig. [7] except that now the different curves are obtained 
by varying v at fixed 7 rather than 7 at fixed v. The 
similarity between Figs. [7J and [TJ] is quite remarkable. 
Depleting the valence level has apparently a very similar 
effect as increasing the bare fermion mass, i.e., it drives 
the NJL2 baryon into the non-relativistic regime where 
it becomes indistinguishable from the GN baryon. 

In order to summarize the results obtained in the 
present work, we have combined all the available points 
from the numerical HF calculations with the analytical 
results taken from Fig. [3] see Fig. [15] The numerical 
points start to fill the white gap in Fig. [3] and match per- 
fectly onto the analytical predictions in the region where 
these can be trusted. This can be seen more clearly in 
the 2d-plots of Figs. Hand Figs-EHU] We do not show 
any results for 7 > 2 where the asymptotic formula (|31j) 
is fully adequate. Notice that when the present work was 
started, only the short solid curve starting at v = 1, 7 = 
had been known from the derivative expansion [l2j. 



chiral limit and to completely filled single particle lev- 
els. After briefly reviewing these results, we have identi- 
fied another parameter region where systematic analyt- 
ical approximations can be performed, namely the non- 
relativistic regime (heavy fermions or weak filling of the 
valence level). Here, a recently developed no-sea effective 
theory has proven to be very efficient. The gap between 
these two asymptotic approaches could only be filled at 
the expense of full numerical HF calculations including 
the Dirac sea. As a result, we have obtained a rather 
comprehensive picture of how the baryon evolves from 
a Skyrme-type topological object in the chiral limit to 
a non-relativistic valence bound state. The transition 
from induced fermion number to valence fermion num- 
ber has a simple interpretation in terms of spectral flow 
as a function of the confinement parameter, as discussed 
in connection with Fig. [6] In the heavy fermion limit, 
it becomes irrelevant whether the model had originally 
a U(l) or Z2 chiral symmetry, and the results for the 
massive NJL2 and GN models converge. 

As compared to the massive GN model with its com- 
plete analytical solution, the analysis of the massive NJL2 
model is significantly more involved. Although the self- 
consistent potentials which we find numerically appear 
to have simple shapes, we have not been able to come 
up with an analytical solution. One source of difficulty 
is the fact that the potentials in the NJL2 model are not 
reflectionless. Related to this, even in those limits where 
we have analytical control, the potentials do not lead to 
any known exactly solvable Dirac equation. If this prob- 
lem is analytically tractable at all, the techniques must 
be quite different from those which have been successful 
in the GN model case. 

Finally, coming back to the phase diagram, we have 
in fact determined one particular phase boundary in the 
present work. Like in the GN model, the T — base line 
of the phase diagram in (7, /i, T) space is expected to be 
given by the curve Mb (7) of Figs. [4] and [5] Another piece 
of information about the phase diagram is the vicinity of 
the tricritical point discussed recently in Ref. [2||. To 
complete this picture remains quite a challenge. 



V. SUMMARY AND CONCLUSIONS 



This paper is part of an ongoing effort to establish the 
phase diagram of (large iV) Gross-Neveu type models. 
Ultimately, we would like to construct the phase diagram 
of the massive NJL2 model in (7, /i, T) space in a similarly 
complete fashion as what has already been achieved for 
the massive GN model. One important building block 
which was still missing are the baryons in the massive 
NJL2 model. The only source of information so far was 
the derivative expansion, restricted to the vicinity of the 
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APPENDIX A: RESULTS FOR COEFFICIENTS 
DEFINED IN SEC. Hill 

Here we collect the 7-dependent coefficients which en- 
ter the results of the no-sea effective theory in Sec. III. 
Coefficients for the scalar potential S, Eq. (|22[) . 



1 



S22 



s 44 



S62 



S64 



«66 



4(1 + 7 ) 2 ' 

3 7 + 4 



S42 



5 7 2 + 3 7 + 2 
967(1 + 7) 5 



647(1 + 7) 4 

91 7 4 - 102 7 3 + 275 7 2 + 276 7 
230407 2 (1 +7) 8 
8I7 3 + 137 2 - 1407 - 24 
9216 7 2 (1 -I-7) 7 
45 7 2 - 136 7 - 60 



(Al) 



9216 7 2 (1 + 7 ) 6 
Coefficients for the pseudoscalar potential P, Eq. (|22|) . 

1 (3 7 + 2)( 7 -3) 



P33 = 



P55 



87(1 + 7) 2 ' 

7~2 
64 7 2 (1 + 7 ) 4 



P53 



192 7 2 (l + 7 ) 5 



(A2) 



Coefficients for the valence fermion density p va i, Eq. 

1 7 2 + 7 + 1 

"22 = ~77~, ! 'I "42 = 



"44 



"62 



"64 



"66 



4(1 + 7)' 48 7 (l + 7 ) 4 

7 + 4 
64 7 (l + 7) 3 

2 7 4 + IO7 3 + 497 2 + 377 + 11 
2880 7 2 (1 + 7 ) 7 
15 7 3 + 1517 2 + 767 + 48 
9216 7 2 (1 + 7 ) 6 
9 7 2 _ 4Q 7 + 12 

9216 7 2 (1 +7) 5 



(A3) 



Coefficients for the induced fermion density p m d, 

Eq. m, 



1 



'42 
&66 



167(1 + 7 ) 3 
2 7 2 -57-4 



?. 44 



ri ? 64 



32 7 (1 + 7 ) 3 
4 7 2 _ 7 7 _ 



1927 2 (l + 7) 6 ' 1287 2 (l + 7) 6 

5(7-2) 



2567 2 (1 +7) 5 



(A4) 
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